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Gate-Drain  Geometry  Effect  on  the  Current-Voltage 

Characteristics  of  GaAs  MESFETS 


1.  INTRODUCTION 

Most  of  tlic  gallium  arsenide  (GaAs)  MESFET  modeling  work  to  date  has 
neglected  to  C)uantify  the  effects  of  gate  to  drain  geometry  on  device 
current-voltage  (I-\')  characteristics*.  A  physical  model  by  Chang  and  Day 
was  recently  reported^  which  includes  the  effects  of  the  gate  to  drain 
distance  I.qq  and  appears  to  agree  well  with  measured  data  for  short  gate- 
drain  spacing  encountered  in  typical  GaAs  FETs.  Recent  interest  in  large 
signal  non-linear  modeling  for  power  amplifier  design  requires  accurate 
knowledge  of  tlie  I-V  characteristics  of  the  MESFETs,  but  relies  on 
empirical  methods  to  obtain  the  I-V  curves^’"*.  An  accurate  computer-based 
model  providing  accurate  DC  characteristic  simulation  based  on  MESFET 
phv’sical  dimensions  and  materials  properties  would  reduce  the  design  time 
of  large  signal  amplifiers.  Although  the  computer  model  presented  here 
uses  measured  data  as  inputs  to  the  simulation  to  explore  the  effect  of 
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*  Pucci.  R.  A.,  Raus.  H.  A.,  and  Statz,  H. 
l>ropcrt,ics  of  gallium  arsenide  microwave  field 
Advances  i  n  El  1  ron  i  cs  and  Electron  Physics.  Vol 
V<'rk.  pp  19o-2(i5. 

^Chang,  C-S  and  Day,  D-Y  S.  (1989)  Analytic  theory  for  current- 
voltage  characteristics  and  field  distribution  of  GaAs  MESFET’s,  IEEE 
Tra us  ■  Elect  ror,  Dev  i  ces  ED-36(No.  2):269-280. 

*llwaug.  V.  1)..  Shih,  V-C..  be,  H.  M.,  and  Itoh,  T.  (1989)  Nonlinear 
modi  1  i  tig  and  verification  of  MMIC  amplifiers  using  the  waveform-balance 
method.  I EEE  Tra  n  s .  M i crowa ve  Theory  and  Techn i nues  MTT-37(No.  12) : 2125- 
2 1  33 . 

■*Curtiee.  V.  II.  (1988)  GaAs  MESFET  modeling  and  nonlinear  CAD.  T EEE 
'  Miorow.o.  I  '■  •  III!  res  VIT  3G(  No .  2  )  :  220  -  230 . 


(1975)  Signal  and  noise 
effect  transistors,  in 
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gate-drain  distance,  the  program  may  be  easily  modified  so  that  physical 
and  material  properties  of  the  GaAs  FET  are  all  that  are  necessary  to 
provitle  accurate  results. 

Although  the  Chang-Day  model  provides  good  results  for  short  Lqq 
distances  up  to  4  micrometers,  we  repoit  a  modified  model  that  is 
experimentally  verified  with  GaAs  FETs  with  L^g  ranging  up  to  12  micrometers 
(pm).  To  determine  the  range  of  validity  of  this  model,  and  to  provide  a 
data  base  of  measured  data  for  testing  improved  models,  this  paper 
describes  the  results  of  a  comprehensive  set  of  GaAs  FET  measurements.  Ten 
MESFETs  with  Lqq  ranging  from  3pm  to  12pm  were  fabricated  and  tested  on  a 
common  semi - i nsulating  GaAs  wafer.  Data  collected  from  DC  FET  parameter 
measurements  were  inputs  to  a  computer  program  based  on  the  Chang-Day  model 
to  plot  simulated  I-V  curves  to  which  measured  I-V  characteristics  were 
compared.  The  model  was  modified  to  incorporate  experimentally  determined 
values  of  pinch-off  voltage  and  to  better  predict  the  I-V  characteristics 
of  longer  Lqq  devices. 

The  Chang-Day  model  was  selected  for  our  MESFET  modeling  because  (1) 
it  incorporates  a  realistic  electron  drift  velocity-electric  field 
relationship®  and  (2)  the  model  simultaneously  solves  the  two-dimensional 
Poisson  equation  and  diffusion  equation  in  closed  form.  Therefore,  a  full 
set  of  I-V  curves  can  be  obtained  within  seconds  based  on  experimental 1\' 
determined  MESFET  parameters.  Another  approach  to  simultaneously  sol\ing 
the  two-dimensional  Poisson’s  equation  and  diffusion  equation  using 
numerical  simulation  was  tried  using  Frensley’s  model.®  Using  this  model 
requires  thousands  of  mesh  points  to  define  the  entire  active  region  of  our 
longer  devices.  These  mesh  points  must  be  constantly  updated  throughout 
tlie  simulation,  leading  to  lengthy  computer  simulation  time.  Numerical 
simulation  clearly  showed  the  high  field  region  near  the  gate  edge  of  the 
ungated  region  between  the  gate  and  drain.  However,  calculated  I-V  results 
vary  widely  from  experimental  data  because  of  the  uncertainties  in  FET 
physical  parameters  such  as  diffusion  coefficients  and  channel  thickness. 

®Chang,  C-S,  and  Fetterman,  H.  R.  (1986)  Electron  drift  velocity 
versus  electric  field  in  GaAs,  Solid  State  Electronics  29: 129.’) -1296. 

®Frensley,  W.  R.  (1981)  Power-limiting  breakdown  effects  in  (iaAs 
MESFET’s,  IEEE  Trans .  Electron  Devices  ED-28(No.  8):962-970. 
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Ve  developed  a  FORTRAN  program  based  on  the  Chang-Day  model  to  derive 
a  full  set  of  I-V  characteristics.  The  program  uses  Chang  and  Day’s 
physically  based  analytical  solution  assuming  a  uniformly  doped  active 
layer.  Results  from  the  computer  model  are  compared  to  a  family  of  MESFETs 
designed  and  tested  at  RADC  and  fabricated  at  a  foundry.  Distance  Lqq  is 
the  only  parameter  varied  among  the  MESFETs,  making  them  particularly  well 
suited  to  modeling  the  I-V  characteristics  as  a  function  of  Lqq,  holding 
all  other  MESFET  physical  parameters  constant. 

2.  DESCRIPTION  OF  EXPERIMENTAL  WORK 

The  family  of  MESFETs  studied  were  fabricated  and  tested  on  a  250pm 
thick  semi-insulating  GaAs  substrate.  The  N-type  channel  is  Si^®  ion- 
implanted  at  a  dual  energy  of  120KeV/30KeV.  The  n"*"  contact  implants  as 
well  as  the  channel  implant  are  activated  by  furnace  annealing  with  As- 
overpressure .  Ni:Ge:Au:Ni  ohmic  metal  was  evaporated  on  the  n"*"  contact 
regions  and  alloyed.  Mesa  isolation  provides  device  isolation.  The 
Ti:Pt:Au  1pm  gate  is  self-aligned  in  a  recessed  channel  etched  for  a 
nominal  pinch-off  voltage  of  2.5  volts.  Polyimide  passivates  exposed 
channel  regions,  and  Ti:PtiAu  overlay  metal  is  evaporated  to  a  thickness  of 
1.2  pm  for  low  resistance  contacts  at  the  probe  pads. 

All  MESFETs  tested  have  6-finger  interdigitated  gate  topology,  1pm 
gate  length  and  600pm  total  gate  width.  Figure  1  shows  the  gate  topology 
and  t.lie  three  segmented  drains  that  evenly  divide  the  gate  periphery. 
Therefore,  the  effective  gate  width  for  each  MESFET  segment  is  200pm.  Each 
MESFET  segment  has  a  unique  gate  to  drain  distance  Lqq  defined  as  the 
distance  between  the  drain-side  edge  of  the  gate  and  the  gate-side  edge  of 
the  drain  as  shown  in  Figure  2.  The  segmented-drain  topology  isolates  Lqq 
as  the  only  physical  parameter  that  changes  wlien  comparing  tlie  I-V 
fdiat  ar  I.er  i  s  (,  i  c'.s  among  MESFET  segments  with  the  same  source,  gate,  Schottky 
bariier,  and  channel  characteristics.  MESFETs  were  designed  with 
ranging  from  3.0pm  to  12.0pm  in  1  pm  increments.  In  addition,  each  MESFET 
contains  one  reference  drain  segment  to  facilitate  matching  devices  with 
similar  pinch-off  voltage,  parasitic  source  and  gate  resistance,  and  built- 
in  voltage?,  again  isolating  the  effect  of  Lqq  on  the  I-V  curves. 
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FIGURE  1.  Topology  of  the  segmented-drain  MESFET  with  three  drains.  Each 
MESFET  segment  has  two  gate  fingers  and  200  ftm  gate  periphery.  The  shaded 
region  designates  source  metal.  Source  metal  is  bridged  over  the  gate. 
Drain  contacts  are  compatible  with  on-wafer  coplanar  waveguide  probe  heads 
with  a  4  mil  pitch.  Note  the  different  gate-drain  distances  among  tlie 
three  MESFET  segments. 
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Figure  2.  (a)  Cross-section  of  MESFET  active  region  defining  the  critical 
dimensions  used  in  the  computer  modeling.  (b)  Electric  field  profile  along 
the  channel.  (c)  Electron  drift  velocity  profile  along  the  channel  sliowing 
field-dependent  drift  velocity  used  in  the  computer  model  (after  Chang  and 
Day,  Ref.  2). 


DC  measurement  techniques  developed  by  Fukui^  were  employed  to 
determine  the  Schottky  barrier  potential  V^,;,  extrinsic  pinch-off  voltage 
Vp,  active  channel  donor  density  Nq,  thickness  a,  and  resistance  Rq?  as  well 
as  parasitic  source,  gate,  and  drain  resistances  IlgM’  ^C. '  ^DM 
respectively.  Expressions  relating  a  MESFET’s  intrinsic  channnel 
characteristics  with  its  extrinsically  measured  parameters  are  essential  to 
modeling  its  I-V  curves.  From  Refs.  7  and  1  the  gate-to-channel  potential 
Wp  required  to  fully  deplete  the  channel  of  carriers  is  given  by 


qN^a^ 
2fs  ’ 


(1) 


where  q  is  the  electron  charge  =  1.602  x  10'^®  C  and  is  the 
dielectric  permittivity  of  the  channel.  For  GaAs,  a  value  of  =  13.1  is 
used  throughout  our  FET  modeling.  The  extrinsically  measured  pinch-off 
voltage  is  related  to  the  channel  pinch-off  potential  by 


Vp  =  Wp  -  Vbi  .  (2) 

Therefore,  the  experimentally  determined  values  of  Vj^j  and  Vp  must  be  added 
before  calculating  the  channel  donor  density  and  thickness.  Although  not 
explicitly  stated  in  the  Chang-Day  model,  Vp  =  Wp  —  Vy  is  incorporated  into 
our  implementation  of  the  model  to  properly  relate  the  terminal  pinch-off 
voltage  with  the  channel  thickness  and  donor  density. 

The  second  equation  required  to  uniquely  determine  Nq  and  a  is  their 
relationship  to  the  measured  saturation  current  Ig: 

Ig  =  qN^aZv^  ,  (3) 

where  Z  if  the  channel  width  =  200pm  for  all  MESFETs  tested  and  is  tlie 
electron  saturation  drift  velocity.  The  saturation  drift  velocity  in 
n— GaAs  is  determined  as  the  temperature  dependent  value® 

^Fukui,  H.  (1979)  Determination  of  the  basic  device  parameters  of  a 
GaAs  MESFET ,  ATfcT  Bell  System  Tecliiiical  Journal  58(No  3);771-T97. 

®Sze,  S.M.  (1981)  Physics  of  Semiconductor  Devices .  2iid  ed . ,  Wiley, 
New  York. 
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wliorc 

and 


VgCT)  =  Vq  -  aT  , 


(4) 


Vq  =  12  X  10®  cm/sec  , 
a  =  0.015  X  10®  cm/sec— K  . 

At  i-ooiu  tempei'ature ,  =  7.5  x  10®  cm/sec,  the  value  used  throughout  our 

MESFET  modeling. 

Determining  Nq  and  a  require  experimentally  determined  values  of 
Vp  and  Ig.  The  Schottky  barrier  built-in  voltage  found  in  Refs.  1  and  6 
is 


V 


bi 


(•5) 


where  A*  is  the  effective  Richardson  constant,  T  is  the  junction 
temperature,  k  is  the  Boltzmann  constant,  and  Jg  is  the  gate  saturation 
current.  Fukui’s  method  is  used  to  find  Jg  by  a  semi-logarithmic  plot  of 
gate  current  as  function  of  gate  forward  bias  and  extrapolating  to  zero 
gate  bias.  Using  A*  =  8.7  Amperes/cm^-K^  to  account  for  the  effective 
electron  mass  in  n-GaAs  at  300  K,  the  equation  for  reduces  to: 


Vy  =  0.350  -  0.0258  InCJg)  . 


(G) 


Tabulated  values  of  determined  for  each  MESFET  are  sliown  in  Table  1. 


Pinch-off  voltage  Vp  is  determined  approximately  b\'  monitoring  the 
drain  source  I-V  characteristics  as  a  function  of  gate-source  bias  \q. 
Drain-source  bias  Vq  is  maintained  at  0.05  volt  to  minimize  depletion 
region  variation  along  tiu'  rlunim-l.  Extrapolating  llie  Ijj  -  plot  to  zero 
drain  current  provides  an  estimate  of  the  pinch-off  voltage  Vp.  This 
initial  value  of  Vp  is  refined  by  a  Fukui  technique  tliat  plots  the 
equivalent  drain-source  resistance  Rpg  =  Vp/Ip  as  a  function  of  parameter  A' 
defined  as 


..-1 


— 

o 

> 

1 

> 

L  ^ 

Vbi  ^  Vpj 
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Tal)le  1:  Measured  and  calculated  MESFET  parameters 


TABLE  1 

DEVICE 

Log  (m-ITi) 

Vbi 

Volts 

Vp 

Volts 

Rsm 

OHMS 

Rdm 

OHMS 

Rg 

OHMS 

Ro 

OHMS 

Nd 

X10’®  Cm^ 

a 

iJLm 

Is 

MA 

FI  445 

3.0 

0.636 

2.39 

2.30 

11.5 

35.5 

5.08 

16.7 

0.1620 

65.0 

FI  445 

4.0 

0.636 

2.39 

2.3 

13.8 

35.5 

5.08 

16.7 

0.1620 

65.0 

FI  555 

5.0 

0.653 

2.28 

1.80 

17.3 

36.4 

5.40 

17.4 

0.1565 

65.2 

FI  665 

6.0 

0.659 

2.45 

m 

20.0 

36.3 

5.50 

15.9 

0.1682 

64.3 

FI  775 

7.0 

0.64 

2.46 

2.85 

24.1 

36.1 

5.34 

16.7 

0.1616 

64.6 

FI  885 

8.0 

0.658 

2.50 

2.33 

27.7 

36.1 

5.29 

16.3 

0.1675 

65.6 

FI  995 

9.0 

0.635 

2.28 

1.44 

30.8 

37.1 

5.71 

16.5 

0.1602 

63.3 

FI  10105 

10.0 

0.64 

2.49 

2.80 

35.6 

36.2 

5.01 

15.7 

0.1699 

64.1 

FI  1115 

11.0 

0.64 

2.49 

1.78 

37.7 

36.5 

5.08 

15.9 

0.1689 

64.5 

FI  12125 

12.0 

■ 

2.30 

1.35 

41.8 

36.7 

5.35 

16.3 

■ 

63.4 

X 

0.646 

2.40 

2.10 

36.2 

5.28 

16.4 

0.1639 

64.5 

a 

0.01 

0.09 

0.51 

0.49 

0.22 

0.50 

0.0044 

0.75 

The  Rqs~^  relationship  is  linear  for  the  correct  value  of  Vp.  The  initial 
value  of  Vp  is  tried  and  adjusted  until  a  linear  relationship  is  found.  The 
linear  Rjjg— JT  plot  yields  the  parasitic  resistances  Rsm»  ^^DM» 

Linear  extrapolation  of  the  plot  to  the  ordinate  yields  Rgj^  +  RpM-  Gate  I-V 
cliaiacterist  ics  are  recorded  as  suggested  by  Fukui  in  Ref.  7  for  a  gate 
current  density  of  lO'*  A/cm^.  For  a  200pm  gate  periphery,  gate  bias 
readings  at  20  mA  gate  current  were  recorded  to  determine  the  effective 
resistance  in  common-source  and  common-drain  configurations.  The 

difference  between  the  common-source  and  common-drain  resistances  yields 
the  quantity  Rg]^  —  RqvI’  Combining  it  with  Rgj^^  +  Rq^j  obtained  from  the 
Rog— X  plot  (Figure  3)  provides  the  measured  values  of  source  and  drain 
resistances  used  in  our  model  and  tabulated  in  Table  1.  The  value  of  Rq;^ 
is  plotted  as  a  function  of  Lqq  in  Figure  4,  showing  a  linear  relationship. 
The  plot  extrapolates  to  the  origin  indicating  a  very  low  ohmic  contact 
resistance.  As  seen  in  Table  1,  Rqm  is  the  distinguishing  parameter  that 
varies  across  the  device  samples. 

The  saturated  channel  current  Ig  is  determined  from  the  drain  I-V 
characteristics  to  be 

Is  =  ip  ~  IsUB  > 


where  Ip  is  the  fully  undepleted  channel  current  evaluated  at  Vg  =  +  Vbi  and 
IgUB  is  the  substrate  leakage  current  evaluated  at  Vq  =  -Vp  with  Vq  at  the 
knee  voltage.  The  value  of  Ig  is  critical  to  determining  Nq  and  a. 
Rearranging  Eq.  (1)  to  solve  for  and  substituting  into  Eq.  (3)  yields 

o 


f'D  =  2r-q-^ 


1 


Vp  + 


I. 


(9) 


with  Vp  =  Vp  +  Vj,;. 


Foi'  n— GaAs,  Eq.  (9)  may  be  simplified  to 

p 

N  -  4.79 

‘'D  -  V  +  Vu- 

'  p  ^  bi 


(10) 


where  is  in  units  of  10*®  cm"*,  Vp  and  V^j  in  volts,  Ig  in  amperes,  and  Z 
in  cm.  In  a  similar  fashion,  the  channel  depth  a  is  found  to  be 


a  =  2(, 


V  W  — 


(11) 


9 


Rds  [OHMS] 


H) 


Rdm  [OHMS] 


Log  (m-hti) 


I’igine  4.  Plot  of  drain  rosistance  as  a  function  of  gatc-ilrain 

distatuu'  showing  a  linear  relationship.  The  plot  is  exti'apolat  «'d  lo 

1  he  or i g i n . 
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For  n— GaAs,  this  expression  simplifies  to 


a  =  0.174(Vp  +  Vbi)  ^  ,  (12^ 

where  a  is  in  microns,  I3  in  amperes,  and  Z  in  cm.  Therefore,  Nq  and  a  are 
derived  from  experimentally  determined  values  of  pinch-off  voltage,  built- 
in  voltage,  and  channel  saturation  current. 

All  experimentally  measured  and  derived  parameters  used  for  our 
MESFET  modeling  are  tabulated  in  Table  1.  Excellent  agreement  in  all 
parameters  is  shown,  with  a  standard  deviation  of  less  than  4  percent  of 
mean  value  for  all  parameters  except  RgM*  large  deviation  of  Rgyj  may 
be  attributed  to  the  graphical  method  of  its  determination.  Rpj^j  is  the 
distinguishing  parameter  among  the  devices  due  to  the  change  in  Lqq  . 

3.  COMPUTER  MODELING 

A  FORTRAN  computer  program  was  developed  following  the  Chang-Day 
model  which  divides  the  channel  into  three  regions  as  shown  in  Figure  2. 
Electric  field  and  electron  transport  properties  are  incorporated  into  the 
solution  of  Poisson’s  equation  and  the  diffusion  equation  in  the  gated  low 
field  region  (region  I),  the  gated  high  field  region  (region  II),  and  the 
ungated  high  field  region  between  the  gate  and  drain  (region  III).  A 
uniform  doping  profile  in  the  active  layer  is  assumed,  vhicli  well 
rcpi'c'sents  the  doping  profile  data  provided  by  the  foundry.  Tlie  uniform 
profile  assumption  simplifies  the  equations  dev'eloped  by  Chang  and  Day  and 
makes  computer  solution  for  tlie  three  regions  particiilarlv'  efficient. 


3.1  Linear  Region 

Following  Eqs.  (42)  and  (43)  in  Ref.  2,  we  calculate  the  channel 
current  for  the  linear  region  by  solving  the  transcendental  equation 


where 

with 


_  q^NpZ/io^ 

6(.L 


^bi  ~  ~ 

W„ 


(I  3a) 
(14) 


and 

( 

ho  \  [Vbi  -  Vg 

+  I^R(^ 

2 

(l.b) 

\ 

“  y  [  '"p 

- 

The 

max i mum 

linear  range  drain  voltage 

''d.l.max 

i  11  Kf| . 

(44) 

of  Ref.  2  is 

il-'cd 

i  n  K(|  . 

(  1  )  )  a  l)o\'e 

and  is  itself  a 

function  of  I 

V.’ 

i  ncrernen  t  1  f. 

uni  i 

1  the  tr 

a  n  scendi'ii  t  a 

1  ef|ua  t  i  on  F|  = 

0  with 

'd  -  'd.kmax 

in  Eq.  (11). 

The 

max imum 

chan  11  el 

current-drain  voltage 

po  i  n  t 

(^C.t.MAX'  '  D.KMAX  ) 

pstalil  i.shes  the  upper  hound  of  the  linear  region. 

Nexi  .  the  interior  I^-Vq  points  of  the  lineal'  region  are  det  eim  i  iied . 
iia  iii<  •  I  y  1  h<  >se  I  ^ j  \  a  I  u  es  that  sa  1  i  s  I  y 

0  <  Ic  <  Ic.LMAX  • 

0  <  Vf)  <  VpLMAX 

I\n  a  gi\.en  \'iilue  of  !(■,  \'q  may  be  obtained  <'xplicill\’  b%’  reuriting  K(|. 

I  Ki  I  i  11  tlie  form : 


as 


function  of  I 


C- 


Since  li,  is  a  fuiicl  ion  of  If.  and  not  Vq,  we  have  a  cubic  ecpiatioii  for 

Solving  for  j  |  explicitly.  we  obtain  a 
corresponding  taliie  of  \'q  by  Eq .  (14).  A  sufficient  number  of  If-Vp)  points 
are  calculated  for  the  I-V  curve  plotting  using  commercial  software 
plotting  routines. 


3.2  Knee  Region 

file  knee  region  is  defined  as  the  cui'ren t -\ o  1  t age  values  satisfying 


^(l.MAX  ^  ^  ^C.KMAX 

'd.i.max  'd  -  ''d.kmax 


(IS) 


with  (Ifi^MAX'  ^OI.M-XX^  determined  fi'om  the  linear  region  solution  as  the 
lowei'  bound.  fhe  next  step  is  to  obtain  the  upper  boiinrl  of  the  knee  region 


'  'cKMAX-  'd.kmax*-  I’ollowing  Ref.  2.  Vp  derived  if  we  let 

,  c,  ... 


'  ::ri 


where 


(a  constant) 


C-’Oa) 


(20b) 


where  Eq  =  vjfi  is  the  ci  itical  electric  field  and  the  maximum  field  in  ilie 
linear  region  is  defined  as 


-  2  [^T  ~  ’ 


where  E-p  is  the  threshold  electric  field  where  the  electron  drift  velocity 
is  maximum.  Rearranging  Eq.  (19),  we  have 


Substituting  into  Eq.  (14)  and  solving  for  Vp  yields  the  expression  for 


'D.KMAX 


''d.kmax  =  '^p  -  Vy  +  Vq  +  Ic.kMAX  *^0  • 

The  other  consequence  of  setting  0i^  =  C'  is  that  U|^  =  0  as  dc.'fined  in  Kc| . 
(49b)  of  Ref.  2.  The  relationship  at 

(48)  and  (46)  in  Ref  2  as 


L  -  (U,  +  UJ  -  [Eq.(46)  of  Ref.  2]  =  =  0 

lnst;i-ting  the  expression  for  Eq.  (46)  of  Ref.  2  and  defining  A, 


,Y  ^  lli!  Is  r 


into  Eq.  (24)  above  will  lead  to  an  expression  for  X  that  may  be  solved  to 
get  IcKMAX-  Observing  that 

r  =  AX  +  B  ,  (201 


( 'Jf-.i  ) 


-  (Vbi  -  Vo) 


(  2lib  ) 
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llicii  from  FC(i .  (14)  with  'd  =  '^D.KMAX^  o"®  obtains 

(  .  (-C'A)  X  +  (1-C'l3)  .  (27) 

Also,  from  I'.ci .  I  IT )  of  Ref.  2 

(  ^  )  =  (-CA)  X  +  (1-CB)  .  (28) 

8ul>- 1  i  t  \i1  in>^  Kci^ .  (27l.  (28),  at\fl  (27))  into  F.q .  (24)  and  performi  np,  some 

ti-dioiis  alpebra  leads  to  the  “pseudo"  eubi<'  ecpialioii 

:\ 

If,  =  by  -t-  b|X  +  b^x'^  +  b^x^  +  l>,,X'  =  0  .  (20) 


Tlie  coefficients  b^^^  are  gi\en  by  long  algebraic  expressions  calculated  in  a 

subroutine  of  the  program.  Equation  (29)  is  solved  using  a  Newton-Raphson 

/  I,  \2 

tout  ini'  with  starling  value  ~  (  "T  )  evalualeil  at  Ic.bMAX'  ^ 
lower  bouiid  of  the  eliannel  current  in  the  knee  region; 


X 


11  +  1 


f2(X„) 


n 


0.  1 . 


(dO) 


It  follows  tliat  tlie  maximum  channel  current  in  the  knee  region  is  gi\en  as 


^C.KMAX 


C^p^finul 


-  V^i  +  Vo] 

Rs 


(31) 


Channel  length  L,^,  where  the  electric  field  is  less  than  E|  .  can  be 
deteimined  by  Eq .  (48)  of  Ref.  2  at  (If-.KMAX’  Vp).;;^l^x).  Channel  length  Lj,, 

where  E,  <  F  <  .  is  calculated  using  Eq.  (52)  of  Ref.  2.  With  1.^^  and 

fiftermined.  t  lie  gated  channel  saturation  length  1,^  is  simpH' 

b,  .  I.  -  1.,  -  1,.  .  (52) 

Willi  the  knee  region  bounds  specified,  the  interior  points  art'  then 
detirmined  by  specifying  a  value  of  Vq  and  calculating  the  cor  resiiontl  i  ng 
\alue  of  I,,  to  '-ol\('  Ft).  (48)  of  Ref.  2  where  ^  0. 


3.3  Saturation  Region 


The  satniation  region  is  defined  as  tlie  current  and  voltagi'  \alues 


"a  t  i s  f y i ng 


>  I  ^ Vv\v 

'd.KMAX  ■  *!)  -  'n,M\V 

1.') 


(53) 


where  V 


D.KMAX 


is  the  terminal  drain  bias  and 


'(  .MAX 


the 


ma  X  1  III  u  m 


attainable  channel  current  for  a  given  (V(^,  bias  condition.  Tliis 

region’s  bounds  are  completely  defined  except  for  IcmaX-  determine  the 

Vq  characteristics  in  the  saturation  region,  we  increment  the  channel 
current  by  AI(^  starting  with  the  lower  bound  Ick.MAX  that 


-  ^C.KMAX  +  AI( 


I  34  ) 


Cliannel  lengths  L^,  and  L2  are  then  calculateri.  TIh-  current  1^-  is 

i  terat.i  vely  incremented  until  becomes  negative  (indicating  a  fii!l>- 

saturated  channel).  ^CMAX  ts  the  maximum  possible  cliannel  current  and 
defines  the  upper  bound  of  the  channel  current  in  the  saturation  region. 

For  each  incremental  value  of  the  ungated  saturation  distance  Lj 

must  be  calculated  to  determine  the  corresponding  V^.  L3  is  obtained  by 
solving  the  transcendental  equation  Eq.  (18)  from  Ref.  2.  Ve  define 
function  fj  such  that 


^3  “  “  [^he  right  hand  side  of  Eq.  (18)  in  Ref.  2=0 


(30) 


Using  a  Newton -Raphson  method,  L3  is  iteratively  determined  using  function 
fg  and  its  derivative  fj'  =  df^/dLj.  A  starting  value  of  I.j  =  10"'*  cm  is  used 
since  it  is  expected  that  the  higli  field  region  length  is  approximately 
this  length®.  With  L3  known,  the  quantity 

A  V  =  V(L3)  -  VX-L^)  (36) 

given  by  Eq.  (19)  of  Ref.  2  can  be  calculated.  It  follows  that 

''(^)  =  '^D.KMAX  +  A  V  .  l'37) 

Careful  investigation  shows  that 

-  '(*'3)  +  IcRp  • 

With  channel  current  and  corresponding  drain  \oltage  \  [j  determined.  Ip 
is  converted  to  Iq  to  account  for  the  substrate  leakage  ciiri  ent  .  Therefore 


III  -  Ip  +  ir- 


(  :!!» ) 


York : p. 344 


®Shur,  M.  .S.  (1987)  CaAs  Devices  and  C  i  ren  its.  Plenum  Press.  New 
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Measured 


wliere  V^/Rp  corresponds  to  the  substrate  leakage  current  IsuB' 

IguB  knee  voltage  yielded  a  typical  value  of  Rp  =  2500  Q,  which 

is  the  value  used  in  this  model. 

The  computer  program  sequentially  finds  the  1-V  relationship  in  all 
three  regions  for  discrete  values  of  Vq.  A  complete  set  of  I-V  curves  can 
be  calculated  and  displayed  within  seconds  with  the  computer  model. 

4.  SIMULATION  RESULTS  AND  MODELING  OF  GATE-DRAIN  REGION 

Figure  5  shows  the  I-V  modeling  compared  with  experimental  results 
foi'  two  MESFLTs  with  radically  different  values  of  (and  thus  Rq)  .  The 

MESFET  with  the  lower  knee  voltage  has  =  dpin.  The  second  FET  has  = 

12/im.  The  values  in  Table  1  specific  to  the  two  devices  are  input 
parameters  in  the  MESFET  model.  Excellent  agreement  between  the  modeled 
and  measured  data  is  evident,  especially  with  the  shorter  device. 

In  order  to  explain  the  effect  that  Lqq  has  on  MESFET  I-V 
characteristics,  we  can  begin  with  the  development  of  Chang  of  Day  in  which 


-  '^DM  +  ’ 

“hore  RnM=  • 

is  the  experimentally  determined  value  of  the  parasitic  drain 
and 


(40) 

(41) 

resi stance 


R(Vg) 


h 

qNpZ/x 


1 

a— h 


1 


(42) 


drain  contact  resistance,  t^  is  the  surface  depletion  width  due 
to  surface  states  on  the  ungated  channel,  t  is  the  total  active  la\er 
thickness,  and  the  depletion  width  h  is 


h 


Wd 


(Vbi  -  Vg) 


(43) 


Notice  that  our  experimental  determination  of  Rqj^^  as  a  futu:tioii  of  E|j(.  is 
accf)nnted  for  in  Eq .  (41).  Figure  4  shows  tliat  R,|^  is  essentially  z<'ro. 

wliich  is  in  agieement  with  drain  contact  data  proxidt'd  b>  tlie  foinidr>. 
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Drain  Current  ID  MA 


F  i  i  5 

with  measured  data.  MESreXs^luh  computer  model  compared 
showing  the  greater  drain-source  bffs  7  """  compared, 
lim)  device.  required  to  saturate  the  longer  (12 
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xo 


Furthermore,  the  slope  of  the  Rdm— Lqq  plot  yields  the  undepleted  channel 
thickness  t  —  t^,  in  the  ungated  region  between  the  gate  and  drain. 
Rearranging  Eq.  (41)  we  obtain 


(t.  -  tj  = 


AL 


DC 


qN^Zp  A  R 


DM 


(4-1) 


Using  the  mean  value  of  Nq  =  1.64  x  10^^  cm'^  and  the  modeling  \alue  of 
/<  =  4.0  X  10^  cm^/v— sec,  we  obtain  a  value  of  (t  —  t^,)  =  0.1375/im.  The 
depletion  width  h  varies  between  0.076pm  at  Vq  =  0  volts  and  0.1530pm  at 
V(^  =  —2.0  volts.  For  the  experimentally  determined  mean  cliannel  thickness 
d  =  0.  1639pm,  R(Vq)  ranges  from  0.15  fl  at  =  0  volts  to  6.1-5  fl  at  \q  = 
—2.0  volts.  It  should  be  noted  that  R(Vq)  rapidly  appi'oaches  very  large 
values  as  Vq  approaches  pinch-off. 

Chang  and  Day  introduce  a  scaling  parameter  5  to  allow  more 
flexibility  in  fitting  the  model  to  measured  data,  in  which  Rq  and  Rg  are 
taken  as 

=  ^DM  + 

•  (46) 


A  \alue  of  ,s  =  1  proved  to  be  adequate  for  our  short  devices,  but  as 
Log  increased,  a  discrepancy  between  the  model  and  measured  data  in  the 
knee  region  developed,  most  notably  for  low  Vq.  Figure  6  clearly  shows 
tills  for  our  =  12pin  device.  The  model  did  not  provide  adequate  channel 
resistance  at  low  gate  bias  to  fit  tlie  data. 

To  remedy  this  problem,  the  channel  resistance  Ro(^'g)  determined 
Ijy  Fukui  in  Ref.  7  is  added  to  tlie  source  resistance,  yielfling 


where 


^sg('g)  - 


-1 


1  - 


Vbi  - 


bi 


(48j 


Wlien  this  additional  term  is  added  to  Rg,  the  model  more  closely  follows  the 
measured  results.  Figure  5  shows  an  overlay  of  the  measured  and  modeled  I- 
V  tdiaracter i s t i cs  for  a  h[)(;  =  4/im  and  a  Eqg  ~  12/im  MESFET.  Tlie  model  shows 
gooil  agreement  with  exp<'r  i  men  1  a  1  data  in  the  linear  an<i  knee  la'gion  for  all 
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Figure  6.  Current-vol tage  simulation  of  a  MESFET  with  LpQ  =  10  fim  showing 
the  deviation  between  the  model  and  measured  data  for  low  gate  bias  in  the 
knee  region.  The  text  describes  a  modification  to  the  model  reducing  this 
deviation  as  shown  in  Figure  5. 
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gate  tiiases.  In  the  saturation  region,  better  agreement  occurs  for  lowt'r 
than  for  higher  .  in  which  the  channel  curi'cnt  is  depressed  due  to 
the  liighci  values  of  lls(;(V(-;)  at  these  gate  biases. 

5.  CONCLUSION 

A  computer  program  based  on  the  physically  derived  Chaiig-Day  model 
was  devf.'loped  to  descrilje  the  drain  ciirren  t-vo  1  t  age  characteristics  f)  f  flaAs 
MlSllils.  Expel  i  mental  analysis  of  a  family  of  ME.SEE'I's  with  various  gate- 
drain  flistances  detfirmined  the  input  parameters  to  tlie  model.  Experimental 
;ind  modeled  results  were  compatu'd  to  si  tidy  t  h<-  model’s  effectiveness  at 
predicting  the  effect  of  varying  gate-drain  distance  on  the  draiti  I-V 
cutves.  The  model  was  modified  to  show  better  agreement  between  measured 
aiul  modeled  result.s  in  the  linear  and  knee  regions  across  the  ftill  range  of 
devices  analy?ced.  The  computer  model  may  be  easilv'  modified  such  tliat  till 
inputs  to  tint  model  are  physical  dimensions  or  quantities,  thereby 
eliminating  the  need  to  perform  Fukui-type  DC  tests  on  devices  before 
simulation.  Designers  of  GaAs  MMICs '  may  then  find  the  model  useful  in 
predicting  the  DC  operating  point  of  integrated  CaAs  MESEETs  for  large- 
signal  design  and  simulation. 
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Rome  Laboratory's  technology  supports  other  AFSC  Product  Divisions,  the 
Air  Force  user  community,  and  other  DOD  and  non-DOD  agencies.  Rome 
Laboratory  maintains  technical  competence  and  research  programs  in  areas 
including,  but  not  limited  to,  communications,  command  and  control,  battle 
management,  intelligence  information  processing,  computational  sciences 
and  software  producibility,  wide  area  surveillance/ sensors,  signal  proces¬ 
sing,  solid  state  sciences,  photonics,  electromagnetic  technology,  super¬ 
conductivity,  and  electronic  reliability/maintainability  and  testability. 


